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Abstract 

Relays in cellular systems are interference limited. The highest end-to-end sum rates are achieved 
when the relays are jointly optimized with the transmit strategy. Unfortunately, interference couples the 
links together making joint optimization challenging. Further, the end-to-end multi-hop performance is 
sensitive to rate mismatch, when some links have a dominant first link while others have a dominant 
second link. This paper proposes an algorithm for designing the linear transmit precoders at the trans- 
mitters and relays of the relay interference broadcast channel, a generic model for relay-based cellular 
systems, to maximize the end-to-end sum-rates. First, the relays are designed to maximize the second-hop 
sum-rates. Next, approximate end-to-end rates that depend on the time-sharing fraction and the second- 
hop rates are used to formulate a sum-utility maximization problem for designing the transmitters. This 
problem is solved by iteratively minimizing the weighted sum of mean square errors. Finally, the norms 
of the transmit precoders at the transmitters are adjusted to eliminate rate mismatch. The proposed 
algorithm allows for distributed implementation and has fast convergence. Numerical results show that 
the proposed algorithm outperforms a reasonable application of single-hop interference management 
strategies separately on two hops. 
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I. Introduction 

Multi-hop communication is one strategy to improve reliability on wireless communication links. The 
idea is to send information through one or more relays, who receive and retransmit from source to desti- 
nation. In cellular system designs that support relaying, the relay is usually a piece of fixed infrastructure 
connected to a power source but not a wired backhaul. Unfortunately, aside from simple repeaters, relays 
have yet to see wide deployment commercial despite a tremendous amount of research [?], [?], [?], 
[?], [?], [?], [?], [?], [?], [?], [?], [?]. A main reason is that relays in cellular systems are sensitive to 
interference [?], [?], [?], [?]. 

In this paper, we propose an algorithm for relaying in multiple-antenna relay interference channels. 
We assume half-duplex decode-and-forward (DF) relays [?] that cannot transmit and receive at the same 
time. Thus, the transmission procedure consists of two stages. In the first stage, the transmitters send data 
to the relays. By assumption, the relays attempt to decode only the signals intended for their associated 
receivers. In the second stage, the relays re-generate and re-encode signals intended to each of their 
associated receivers before retransmitting to the receivers. While each receiver intends to receive data 
from only a single transmitter via a single relay, each transmitter can have independent data for multiple 
receivers. Consequently, each transmitter may ask for the aid of multiple relays at the same time and 
each relay may simultaneously forward independent data from a transmitter to multiple receivers. 

Using shared radio resources, the transmissions from the transmitters to relays interfere with each 
other. Similarly, the transmissions from the relays to receiver interfere with each other. If considered 
separately, each stage of transmission is an instance of the single-hop interference broadcast channel [?], 
[?], a generalization of the conventional single-hop interference channel. Each transmitter in the broadcast 
channel has independent data for multiple receivers while each transmitter in the conventional channel has 
data for only one receiver. Although recent results on the single-hop interference broadcast channel [?], 
[?] can be applied separately for the two hops, in this paper we show that higher end-to-end sum rates 
can be achieved when the transmitters and relays are configured jointly. Unfortunately, jointly configuring 
the transmitters and relays is challenging, especially with limited information about the interferers. 

In general, the optimal transmit and receive strategies for sum-rate maximization in interference 
channels are not widely known, even for single-hop channels. Thus, we adopt a pragmatic approach 
that treats interference as noise and maximizes end-to-end sum-rates by searching within the class of 
linear transmit and receive strategies. Assuming the receivers always use the optimal linear MMSE receive 
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filters, we focus on designing the transmit precoders at the transmitters and relays. Unfortunately, the 
problem is nonconvex and NP-hard. Moreover, by definition, the end-to-end achievable rates of two-hop 
links are not continuously differentiable at every point. Thus, finding the stationary points of the problem, 
including its globally and locally optimal solutions, is challenging [?], [?], [?]. 

We assume all two-hop links have a common timesharing value, i.e., the same fractions of time 
for transmission on two hops because decode-and-forward is assumed. The achievable end-to-end rate 
corresponding to a relay is defined as the minimum of the achievable normalized rate from its associated 
transmitter to itself and the achievable normalized sum-rates from itself to its associated receivers [?]. A 
two-hop rate mismatch occurs when some links have a dominant first hop while others have a dominant 
second hop, resulting in low end-to-end sum-rates. An efficient system design should not cause any 
two-hop rate mismatch while mitigating interference. 

Transmit precoder design has been studied widely for the multiple-antenna single-hop interference 
broadcast channel [?], [?], especially for its special case of the single-hop interference channel [?], 
[?], [?], [?], [?], [?], [?]. The methods in [?], [?], [?], [?] were based on the so-called interference 
pricing framework where the transmitters configure themselves based on interference prices fed back 
from the receivers. Interference prices represent the marginal decrease in the sum-utility function per unit 
increase in interference power. In [?], [?], [?], based on a relationship between mutual information and 
mean squared error (MSE), sum-utility maximization problems were solved via iterative minimization 
of weighted sum-MSEs. Under certain conditions on the utility functions, the algorithms in [?], [?], 
[?], [?], [?], [?], [?] are guaranteed to converge to the stationary points of the corresponding sum- 
utility maximization problems. Note that the existing single-hop results are designed specifically for the 
single-hop interference channel. It is not straightforward to incorporate the special features of the relay 
interference channel, e.g., relay signal processing operation and multi-hop transmission. 

To the best of our knowledge, there has been little prior work on transmit precoder design in the 
relay interference broadcast channel. Prior work focused on interference mitigation for special cases of 
this DF relay model [?], [?], [?], [?]. In [?], the authors proposed a transmit precoder design for the 
MIMO relay broadcast channel where a single MIMO DF relay forwards data from a single transmitter 
to multiple receivers. This means that there is no interference on the first hop and there is no inter-relay 
interference on the second hop. Prior work in [?], [?], [?] considered the DF relay interference channel 
where the receivers are equipped with a single antenna and each relay is dedicated to aiding a single 
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transmitter-receiver pair. The transmit precoders at the transmitters and relays are designed completely 
independently in [?] to maximize the minimum average signal-to-interference-plus-noise ratio (SINR) on 
each hop. Prior work in [?] investigated the end-to-end sum-rate performance of different relay cooperation 
strategies and found that the zero-forcing based approach is attractive for interference management. Based 
on the interference pricing framework, the algorithm in our prior work in [?] used approximations of 
end-to-end rates to compute interference prices for designing the second-hop transmit precoders with 
fixed first-hop transmit precoders. It is not straightforward to extend the results of [?] to the general 
relay interference broadcast channel with multiple-antenna receivers. There have been other algorithms 
for designing transmit precoders at the relays and/or transmitters in the relay interference channel [?], [?], 
[?], [?], [?], [?], [?], [?], [?], [?], [?], [?], [?], [?], [?]. Nevertheless, much of the prior work considered 
either AF relays [?], [?], [?], [?], [?], [?], [?], [?], [?], [?], [?], [?] or other relay architectures, like 
the shared relay [?], [?] or two-way relay [?]. In addition, in this prior work, each relay simultaneously 
forwards data for multiple transmitter-receiver pairs unlike in our approach. 

In this paper, we propose a cooperative algorithm for efficiently finding suboptimal solutions of 
the transmit precoder design problem with high end-to-end sum-rates. The proposed algorithm can be 
implemented in a distributed fashion with low communication overhead. The algorithm consists of three 
phases in the following order: i) second-hop transmit precoder design, ii) first-hop transmit precoder 
design, and iii) first-hop transmit power control. In the first phase, we ignore the first hop and focus on 
configuring the relays to maximize the achievable second-hop sum-rates. Essentially, the second hop is 
treated as the conventional single-hop interference broadcast channel. Thus, existing single-hop algorithms 
can be applied to find the stationary points of second-hop sum-rate maximization [?], [?], [?], [?], [?], 
[?], [?], [?], [?], [?], [?], [?], [?], [?]. Having computed the second-hop transmit precoders, each relay 
computes the sum of achievable rates from itself to its associated receivers, which is used as input to the 
second phase. 

The second phase focuses on designing the first-hop transmit precoders. In the naive approach, this can 
be done by applying the prior work in [?], [?], [?], [?], [?], [?], [?], [?], [?], [?], [?], [?], [?], [?] while 
ignoring the designed second-hop transmit precoders. The naive approach, however, may cause a two-hop 
rate mismatch because it does not take into account the timesharing value and second-hop configuration. 
To overcome this limitation, we formulate and solve a new problem to maximize an approximation of 
the achievable end-to-end sum-rates. Such approximations of the achievable end-to-end rates depend not 
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only on the first-hop transmit precoders, but also on the timesharing value and second-hop configuration. 
This allows for second-hop interference mitigation at the same time as rate-mismatch alleviation. Some 
guidelines for selecting such approximations are provided. Having defined a more comprehensive utility 
function, we use the technique in [?], [?], [?] to develop an iterative method that is guaranteed to converge 
to the stationary points of the new sum-utility maximization problem. This concludes the second phase 
of the proposed algorithm. 

The output of the second phase may contain a residual two-hop rate mismatch since only an approximate 
solution is proposed. In the final phase, we propose to fix the shapes of the first-hop transmit precoders 
and to adjust their norms to eliminate completely two-hop rate mismatching. Essentially, this is a transmit 
power control problem. Note that for a two-hop link with a dominant first hop, excess power is allocated 
for the transmissions on the first hop. We propose a method for simultaneously reducing excess power 
for the first-hop transmissions so that the achievable end-to-end rates for all the relays are nondecreasing 
over iterations, thus potentially improving the achievable end-to-end sum-rates. The method is guaranteed 
to converge. At the convergence point, there are no two-hop links with a dominant first hop, i.e., there 
is no longer two-hop rate mismatch. Although there have been many power control algorithms for the 
single-hop interference channel [?], [?], [?], [?], [?], [?], [?], [?], [?], [?] and for the relay interference 
channel [?], [?], [?], [?], they are not designed to eliminate two-hop rate mismatching. Therefore, even if 
applicable to the third phase, existing power control algorithms may worsen the two-hop rate mismatch 
situation. 

We use Monte Carlo simulations to evaluate the average achievable end-to-end sum-rates of the 
proposed algorithm for various relay interference broadcast channel configurations. The naive approach 
of applying existing single-hop results separately for the two hops is selected as the baseline strategy. 
The proposed algorithm and the naive approach have the same second-hop transmit precoders. While the 
timesharing value and second-hop configuration are taken into account in the first-hop transmit precoder 
design in the last two phases of the proposed algorithm, they are ignored in that of the naive approach. 
Numerical results show that the first two phases of the proposed algorithm are enough to provide large 
end-to-end sum-rate gains over the naive approach. In addition, the last phase of the proposed algorithm 
makes considerable improvements in end-to-end sum-rate performance over the output of the second 
phase. This highlights the importance of two-hop rate matching in the DF relay interference (broadcast) 
channel. Finally, each phase of the proposed algorithm is observed to converge in a few iterations. Note 
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that the proposed algorithm can be implemented in a distributed manner with a little more overhead than 
the naive approach. Thus, it is suitable for practical implementation in DF relay networks. 

The organization of the remainder of this paper is as follows. Section In] describes the system model. 



Section III formulates the design problem and discusses the challenges. Section IV presents the proposed 



algorithm in detail. Section [V] numerically evaluates the achievable end-to-end sum-rates of the proposed 



algorithm. Section VI concludes this paper and suggests future research. 

Notation: We use normal letters (e.g., a) for scalars, lowercase and uppercase boldface letters (e.g., h 
and H) for column vectors and matrices. Ijy is the identity matrix of size N x N. For a matrix A, A* 
is the conjugate transpose, tr(A) the trace, |A| the determinant, and ||A||f the Frobenius norm. E[-] is 
the statistical expectation operator. ()( n ) denotes iteration index. ()x is used for transmitters' parameters, 
()r for receivers', and Qx for relays'. 

II. System Model 

Consider a relay interference broadcast channel where K^ transmitters communicate with Kr receivers 
with the aid of Kx half-duplex DF relays, as illustrated in Fig. [T] Each transmitter is assigned a unique 
index from /Ct — {1, • • • , -Kt}- Similarly, each relay is assigned a unique index from /Cx = {1, • • • , Kx} 
and each receiver is assigned a unique index from /Cr — {1, • • ■ , K-r}. Each transmitter may require 
the aid of multiple relays to simultaneously send independent data streams to its receivers. Each relay is 
dedicated to serving multiple receivers that are associated with a single transmitter. Each receiver intends 
to receive data from only one transmitter with the aid of a single relay. Let x(fc) E /Cx denote the index 
of the relay that aids receiver k E /Cr. Let ^i(k) E /Ct denote the index of the transmitter that is aided 
by relay k E /Cx- The transmitters and relays do not share data. We assume that each relay k does not 
attempt to decode the signal intended for receiver m E /Cr with x(m) ^ k. Transmitter k E /Ct has 
A^T,jfc antennas, relay m E /Cx has Nx, m antennas, and receiver q E /Cr has N^ q antennas. 

Half-duplex relays cannot transmit and receive at the same time, thus the transmission procedure 
requires two stages. Using a common frequency at the same time, the transmissions in the same stage 
interfere with each other. For tractable analysis, we assume Gaussian signaling is used in both stages 
although it may not be optimal for the relay interference broadcast channel. In the first stage, the 
transmitters send data to the relays. Treating unwanted signals as additive Gaussian noise, each relay 
decodes its desired signal. Each relay separates its decoded signals for individual associated receivers, 
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Fig. 1. A relay interference broadcast channel where a number of half-duplex DF relays aid the data transmission from a 
number of transmitters (TXs) to their associated receivers (RXs). The solid lines connect the associated nodes and represent 
communication links. The dashed lines represent interference links. 



re-encodes and retransmits to its associated receivers in the second stage. Each receiver treats unwanted 
signals as additive Gaussian noise when decoding its desired signal. 

We consider slowly- varying, frequency-flat, block-fading channels. We denote H m ^ G C x '™ x T >' t as 
the matrix channel between transmitter k and relay m for k G /Ct and m G /Cx- Let Xi^ G C dl fcXl denote 
the symbol vector that transmitter fi(k) G /Ct sends to relay k G JCx in the first stage, where is the 
number of data streams and E(xi ^x^ . ) = 1^ k . We denote F^jt G C NT ^ (k)Xdl - k as the linear transmit 
precoder that transmitter fi(k) uses to send xi to relay k G /Cx- We define Ft = {Ft i, • • • , Ft,k x } G 
Ft = C iVT '"W xdl ' 1 x • • • x C^i'xi^^x. Let p Ttk be the sum transmit power at transmitter k G /C T - 
The sum transmit power constraint at transmitter k G /Ct is 

J2 H F T,m|||= J] tr (F T , m F Tjm ) <PT,fc,VfcG/C T . (1) 
fj,(m)=k /j,(m)=k 

We denote nx,fc G C Arx fcXl as spatially white, additive Gaussian noise at relay k G /Cx with E(nx,A;n x fc ) = 
°x fc^x Relay G /Cx observes 

yx.fc = H fciAt ( fc )F T ,fcXi ifc + ^ H feiAt ( m )F T , m xi im +n x ,fc. (2) 

v " . r~ — ~^ meCx 
desired signal mj^k 



interference 
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The interference plus noise covariance matrix at relay k G /Cx is 



Rx,fc(F T ) = ^2 H fc,MM F T,m F T,m H fc, M (m) + cr X,fc I Afx, fc - ( 3 ) 
m6/Cx 

Each relay k G /Cx applies a linear receive filter W x ,k G C JVx -* xdl '* to y Xfe . We define W x = 
(W x ,i, • • • , W x ,if x ) G w x = C JVx - lXdl ' 1 x • • • x C^^xX^.^x. The maximum achievable rate on the 
first hop from transmitter fi(k) G /Ct to relay G /Cx is 

i? l!fe (F T ) = log 2 det (i + F* Tik U* kMk) [K^FtT 1 H fejM(fe) F T , fc ) , (4) 

which can be achieved using the linear MMSE receive filter 

1 _1 

H fc,/u(*0 F T,fc F T,fc H fc, |ll (fc) + Rx,A;(Ft) H fejjU ( fe) F T)fc . (5) 



W MMSE 



We define £i jfc (F T ) = 2 Rl - k ^ - 1. Intuitively, £i jfc (F T ) can be thought of as the effective first-hop 
SINR to relay k if we send a single data stream from transmitter n(k) G /Ct to relay k at the rate 

fli,fc(F T ). 

Let G m fe G C iYR " ,xArx fc denote the matrix channel between relay fc G /Cx and receiver m G /Cr. Let 
X2, m £ x<m)Xd2 m denote the transmit symbol vector that relay x{ m ) G sends to receiver m G /Cr 
where cfo.m is the number of data streams and E(x 2;m X2 m ) = Lj 2m . We denote Fx,™ G C JVx '* (m)Xd2 ' m 
as the linear transmit precoder that relay ~x,(rn) uses to send ~x-2,m 

to receiver m G /Cr. We define 

F x = {F x ,i, • • • ,F x ,tf R } G F x = C JV *'*< 1 > xd2 ' 1 x • • • x C^*-***' Let p x ,fc be the sum transmit 
power at relay k G /Cx- The sum transmit power constraint at relay k G /Cx is 

]T ||Fx, m |||= £ tr(F x , m F Xi J <px, fc ,VfeG/C x . (6) 

X(m)=k x(m)=k 

We denote nR m G C iYR < mXl as spatially white, additive Gaussian noise at receiver m G /Cr with 
E(n Rim n Rm ) = <j\ m l NK m . Receiver m G /Cr observes 

VRm = G mjX ( m )Fx,mX2, m + ^ G m ^ q )Fx,q^-2,q +nR >m . (7) 

desired signal «g£ 



interference 
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The interference plus noise covariance matrix at receiver m G /Cr is 

R-R,m(F X ) = ^2 G m)X ( 9 )Fx, ? F X)(? G^ )X(9) +Cr| !m lAT R , m . (8) 

Each receiver to G /Cr applies a linear receive filter WR m G C^'™^ 2 '"* to yR im . The maximum 
achievable rate at receiver to G /Cr is 

^2, m (F x ) = log 2 det (i + F Xjm G^ iX(m) ^^(Fx)]- 1 G mjX(m) F x , m ) , (9) 
which can be achieved by the following linear MMSE receive filter 

-l 

G miX ( m )F x , m . (10) 
The sum of maximum achievable second-hop rates for relay k G )Cx is defined as 

«2, fc ,sum(Fx) = ^, m (F X ). (11) 

X(m)=k 

We define 6,fc(F x ) = 2 i?2 fc — ( Fx ) - 1. Similar to £i, fe (F T ), we refer to f 2) fc(F x ) as the effective SINR 
on the second hop corresponding to relay k. 

We assume the transmissions in each stage start and end at the same time. Let t be the fraction of time 
for transmission on the first hop, which is also referred to as the timesharing value. The fraction of time 
for transmission on the second hop is (1 — t). For example, in 3GPP LTE- Advanced cellular systems, t 
depends on the number of subframes allocated to the backhaul links (i.e., between base stations and relays) 
in a radio frame [?]. In this paper, t is a fixed parameter. The optimization of the timesharing value is left 
for future work. For relay k G )Cx, the normalized rate on the first hop is ^i ^Ft) while the normalized 
sum of second-hop rates for this relay is (1 — t)-R2,fc,sum(Fx)- Based on the relative comparison of the 
normalized rates on two hops, the two-hop link corresponding to relay k G /Cx can be classified into 
the following three categories: i) first-hop dominant if ^^(Ft) > (1 — t)^2,fc,sum(Fx), ii) second-hop 
dominant if tR 1:k (F T ) < (1 - t)R 2 ,k,snm(Fx), and iii) equal rate if tR ltk (F T ) = (1 - t) J R 2 ,fc,sum(Fx). 

The achievable end-to-end rate for relay k G )Cx is defined as the minimum of the normalized 



W 



MMSE 
R.m 



G m ,x(rra)Fx. 



m F X,mGm ]X ( m ) + R-Rm( 
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achievable rates on two hops [?] 

R k (F T ,F x ) = mm{tR hk (F T ),(l-t)R 2Asnm (F x )} (12) 
= min{tlog 2 (l + a,fe(F T )),(l-t)log 2 (l + 6 ) it(FR))}. (13) 

Given the achievable rates on two hops, the end-to-end achievable rate for relay k can be obtained by 
properly allocating fractions of the normalized first-hop rate ti^ ^Fx) for sending data for its associated 
receivers. Specifically, let f3 k . m > be the first-hop rate corresponding to receiver m £ /Cr such that 
x(m) = k, where 

Pk,m<tRl, k (F T ). (14) 

m£K. R 

X(m)=k 



For example, we can achieve R k (FT,Fx) in (12i if relay k uses the following first-hop rate allocation 
strategy 



(1 - t)R 2>m (Fx), if tR ltk (F T ) > (1 - t)R 2tktSum (F x ) 

(15) 

^ 1)fc (F T ) ^gfe , otherwise. 



We can check that fi k ^ m for k £ /Cx and m G /Cr satisfies the condition in ( 14 1. The end-to-end sum-rate 
is defined as 

(Ft, Fx) = t,Fx). (16) 

fce/Cx 

The design of Ft and Fx for maximizing i? sum (FT,Fx) should take into account t. 



III. Problem Formulation 

This section formulates the design problem and discusses the challenges in finding optimal solutions. 
The problem of designing the transmit precoders at the transmitters and relays to maximize the sum of 
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achievable end-to-end rates is formulated as follows 



(OV) : 



max 

(F T ,Fx)GF T xFx 



-Rsum(FT, Fx) 



subject to 



£ tr(F x>m F Xim ) < p x>k , Vk € JC X 



(17) 



x (m)=k 



£ tr(F T)m F^ m ) < p T , fc ,Vfc e /C T . 



(18) 



fi(m)=k 



The transmit precoder design problem for sum-rate maximization in the single-hop interference channel is 
nonconvex and NP-hard [?], [?]. This means that its globally optimal solutions cannot be found efficiently 
in terms of computational complexity even in a centralized fashion. The more complicated problem, (OV) 
is expected to be NP-hard as well. Moreover, (OV) is a sum-utility maximization problem where the 
per-user utility function is not smooth at every point. Thus, it is challenging to find the stationary points 
of (OV), including its optimal solutions [?]. 

In this paper, we focus on finding suboptimal solutions to (OV) with high values of achievable end- 
to-end sum-rates. Remark [1] and Remark [2] discusses two main challenges in solving for high-quality 
suboptimal solutions of (OV). 

Remark 1: Interference mitigation is a challenge in end-to-end sum-rate maximization. According 
to ([2]), each relay observes undesired signals from unintended transmitters on the first hop. Similarly, 
according to Q, each receiver observes undesired signals from unintended relays as well as from its 
associated relay (but they are intended for other receivers). Due to interference, there exists coupling 
among the achievable rates on the same hop. 

Remark 2: Two-hop rate matching is a challenge in maximizing the end-to-end sum-rates of the DF 
relay interference broadcast channel. Specifically, for a given t, Ft, and Fr, there may exist a mismatch 
between the normalized achievable rates on two hops. By definition, a two-hop rate mismatch occurs 
when there exist k,m£ Kx and k ^ m, such that tR\ fe(Fx) > (1 — i)-R2,fc,sum(Fx) and tRi !Tn (F^) < 
(1 — i)i?2,m,sum(Fx)- When this happens, it is always possible to improve the end-to-end sum-rate 
performance of the system design. For example, the norm of Fj & can be scaled down to obtain a 
new set of transmit precoders (F T ,Fx) so that tRi t k(FL) = (1 — t)i?2,fc,sum(Fx)- This decreases the 
interference power from transmitter k to all other relays on the first hop, improving the achievable rates to 
all other relays, especially ti?i i7n (F T ) > tRi >m (F^). This means that i? sum (F T ,Fx) > -R S um(FT, Fx). 
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Thus, an efficient transmit precoder design in terms of end-to-end sum-rate maximization should not 
cause any two-hop rate mismatch. 

IV. Transmit Beamforming Design 



This section presents an algorithm for finding high-quality suboptimal solutions to (OV). Section IV-A 



discusses briefly the design of Fx- Section IV-B presents the design of Ft- 



A. Second-Hop Transmit Beamforming Design 

We design Fx by treating the transmission on the second hop as the single-hop interference broadcast 
channel. Specifically, we need to solve the following optimization problem 

(SV) : max V] -R2,fc,sum(Fx) 

subject to tr (Fx, m F x , m ) < Px,k, V& G JC X - (19) 

x (m)=k 

Note that (SV) is nonconvex and NP-hard. Nevertheless, its stationary points can be found by existing 
algorithms for the single-hop interference broadcast channel. The principle of many existing algorithms 
is to formulate a series of related optimization problems that can be solved in polynomial time by 
available methods and provide multiple approximations or relaxations of the original sum-utility problem. 
In general, the globally optimal solutions of these related problems converge to the stationary points of 
the original sum-utility maximization problem. The key requirement for the applicability of existing 
algorithms is that the utility function of the original problem is continuously differentiable at every point. 
An example is the algorithm for transmit precoder design via matrix-weighted sum-MSE minimization 
in [?], [?]. Due to space limits, the details of the algorithm are deferred to [?], [?]. 

Let Fx denote the resulting second-hop transmit precoders. We denote £ 2 ,k = £2 fc(Fx) for k G ICx- 
By setting i? 2 ,jfc,sum(Fx) equal to log 2 (l + 6,fc), we obtain 

6,* = 2^.*.—^) - 1. (20) 

We assume that each receiver m G /Cr feeds back the value of i? 2 , m (Fx) to its associated relay, i.e., 
relay x( m )- Then, we assume that each relay k G Kx can compute i?2,fc,sum(Fx) from its corresponding 
-R2,m(Fx) such that k = x( m ) an d then &,k- 
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B. First-Hop Transmit Beamforming Design 

1 ) Subproblem Formulation and Challenges: This section focuses on designing Ft given knowledge 
of t and £ 2 ,fc for k G /Cx- It follows from (OV) that the problem for designing Ft is 

(TV) : max V min{i log 2 (l + £i, fc (F T )), (1 - t) log 2 (l + £ 2 , fc )} 
f t sf t z — ' 

subject to tr ( F T,mF Tim ) < pr,k, G /C T . (21) 

/i(m)=fe 

Note that (TV) belongs to the same class of nonconvex and NP-hard sum-utility maximization problems 
as (SV) and it is even more complicated than (SV). While the objective function of (SV) depends only 
on the corresponding transmit precoders, i.e., Fx, that of (TV) depends not only on Ft but also on t 
and ^ for k G /Cx- 

Remark 3: It is not possible to apply existing algorithms developed for the single-hop interference 



broadcast channel to find the stationary points of (TV). As discussed in Section IV-A existing algorithms 
require that the utility function of sum-utility maximization problems be continuously differentiable at 
every point. Due to the min operation, however, the utility function of (TV) is not continuously differen- 
tiable with respect to £i.a;(Ft) at me point that makes t log 2 (l + £i,/c(Ft)) equal to (1 — t) log 2 (l + &,k)- 

Remark 4: In the naive approach, the timesharing and second-hop configuration are ignored, leading 
to an approximation of the objective function Y^keK 1°§2(1 +£i,&(Ft))- The resulting problem has the 
same form as (SV). Thus, its stationary points can be found by existing single-hop algorithms. 

For notational convenience, we define r/k as the following function of &,k and t 

r}k = (l+hk)^ "I- (22) 

Note that tlog 2 (l + rj^) = (1 — t)log 2 (l + &,k)- This means that Ci^Ft) is equal to rj^ when the 
achievable first-hop rate at relay k matches with the sum of achievable second-hop rates from relay k to 
its associated receivers. Thus, rjk is the rate-matching received SINR at relay k. 

2) Proposed Approach: It is challenging to find the stationary points of (TV) because the utility 
function in the problem is not continuously differentiable at every point. In the section, we aim at finding 
suboptimal solutions to (TV) with high end-to-end sum-rates. Instead of solving directly (TV), we 
propose to formulate and solve a new sum-utility maximization problem, which we refer to as (ATV). 
Having the same constraints as (TV), (ATV) uses an approximation of min{tlog 2 (l + ^i fc(FT)), (1 — 
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i)log 2 (l + &,k)} as the utility function. Note that such approximations depend on not only Ft but 
also t and &,k for k S /Cx- Let (Ft), t, ^2,k) denote the utility function of (AFV). In addition, 

we propose to solve (AFV) via iterative minimization of weighted sum-MSEs, the well-established 
technique that has been used widely in prior work [?], [?], [?], [?]. 

Some guidelines for selecting u k (^x )k (F^), t, ^2,fe) are provided. First, it must be twice continuously 
differentiable with respect to £i,fc(ITT) at every point. Second, it must satisfy the following condition 

(l + ei,fc(FT)K(Ci,fc(FT),t,e2,fc) + 2«j fc (ei ) fc(F T ) > t,&,fc) > 0. (23) 

This condition is required to ensure that the resulting iterative algorithm for solving (AFV) is convergent 
as shown in Proposition 1 in [?]. There are many approximate functions that satisfy the conditions in 
the guidelines. It is still unclear, however, what is the best approximation. 

According to Remark [5J the utility function of (FV) is not continuously differentiable with respect to 
£i,A:(Ft) at Ci,/c(Ft) = f]k- We propose an approximation of the utility function of (FV) as follows 

tlog 2 (l + £i,fc(F T )), if 6,fc(F T ) < Vk, 
Uk(£l,fc(F T ), = { (24) 

t log 2 (l + rj k ) + ^3 exp (l - T+g^y) - 1 , otherwise. 

Note that if £i&(Ft) < %, then u k (£i fc(Fx), t, £ 2 ft) is exactly equal to tRi^F^) and hence equal to 
i4(F T ,F x ). If Ci,fe(F T ) > Vk, then u k (^ k (F T ),t,^ 2 ,k) is larger than both tR hk (F T ) and R k (F T ,F x ). 
Moreover, when £i,&(Ft) > rj k , the gap between ^^^(Ft). I2.it) and i?fc(FT,Fx) increases with 



£i,A:(Ft) and is upper bounded by i/log2. Using the approximation Ufc(£i,ft(FT), t, ^2,k) in (24 1 as the 
utility function, we formulate (AFV) as follows 

(AFV): max V «*Ki,fc(F T ), 

subject to ^ tr(F Tjm F Tim ) < p T ,ft, Vfc e /C T . (25) 

/x(m)=fc 

The stationary points of (AFV) are expected to correspond to suboptimal solutions to (FV). 

3) Sum-Utility Maximization via Matrix-Weighted Sum-MSE Minimization: We develop an algorithm 
for solving (AFV) via an iterative minimization of weighted sum-MSEs. We denote E^Ft, Wxfc) as 
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the MSE matrix at relay k G /Cx, which is given by 



E fe (F T ,W Xj( 



H fc,^) FT .fe F T,fc H fc,M(fe) +H-x,fc(F T ' 



W x „ 



tr(W Xifc H fciM(fc) F T , fc ) - tr(F^ fe H* >Mfc) W X)fe ) + d ljJfc 



The MSE of the estimate of Xi^ based on W x fc yx,fc is given by 



e fc (F T ,W x „ 



E||x ljfe - W Xife y x ,fc||| 
tr(E fe (F T ,W x , fc )). 



(26) 



Fixing Ft and solving l9Efc g^'» Wx, ' c ^ = 0, we can check that e k (F^, Wx,fc) is minimized by the linear 
MMSE receive filter W^ 13 given in Q. We denote E M (F T ) = E fc (F T , W™ SE ). After some 
manipulation, we obtain the following well-known relationship between ^^(Ft) and E^ (Ft) 



i2i,*(F T ) =-log 2 (det(E M (F T ))), 



(27) 



or equivalently 



ei, fc (F T ) = l/det(E fci o(F T ))-l. 



(28) 



Based on (28 1, we define 



g k (B kfi (F T )) 



log 2 



u k (l/det(E M (F T ))-l,t,e2„ 



log(det(E M (F T ))), if det(E M (F T )) > (1 + ?7fc) , 

- log(l + r/ k ) - exp [!-(! + Vk) det(E fe (F T ))] + 1, otherwise. 



(29) 



(30) 



We can check that <7fc(Efc,o(FT)) is twice continuously differentiable with respect to Efc i o(Fx) at any 
point. Moreover, g k (E k (Ft)) is a strictly concave function of Efe j o(Fx). Note that (AFP) is equivalent 
to the following optimization problem 

(grP) : nun £ g k (E kfi (F T )) 



fc&Kx 



subject to ^ tr(F Tjm F Tim ) < p T)A , Vfc £ /C T . 



(31) 



me/Cx 
fj,(m)=k 
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Note that (QFV) is nonconvex and NP-hard. We define the following function 

1, if x > (1 + rjk)- 1 

a(x) = < 

(1 + Vk)x exp [(1 + rj k )x — 1] , otherwise. 
The first-order gradient of g k (E k fl(F^)) with respect to (Ft) is given by 

Vg k (E kfi (F T )) 



(32) 



a dg k {E kfi {F T )) 



<9E fcj0 (F T ) 

= a(det(E fci0 (F T )))(E fc ,o(FT))" 1 . (33) 

According to Theorem 2 in [?], the inverse mapping of Vg k (F kt o(FT)) is well-defined. We refer to it 
as 7 fc (-) : <C dl ' fcXdl ' fc C dl - kXdl - k . 

We now use the technique in [?], [?] to solve {AFT) via matrix-weighted sum-MSE minimization. 
We introduce auxiliary variables V fc G C dl - kXdl - k for k G /C x . We define V = (Vi, • • • , V Kx ) G V = 



X • • • X 



(j~<di iKx xdi,A' x We define the following matrix-weighted sum-MSE function 



Sfe (F T ,W x ,V) = tr(Vj£E fc (FT,Wx, fc )) + 9k(7k(Vk)) - tr(V^) 7fc (V fe ). (34) 
A matrix-weighted sum-MSE minimization problem is formulated as follows 
(MTV): min Vs fc (F T ,W x ,V) 

(F T ,W X ,V SFtXWxXV /—i 



subject to ^2 tr( F T, m F Tm ) < PT,fc,V& G ^t- 



/i(m)=fc 

It follows from Theorem 2 in [?] that (MFV) and (QFV) have the same stationary points if the relays 
use their corresponding linear MMSE receive filters and the matrix weights V k = Vg k (E k (FT, W x )) 
for any Ft and W x . Thus, we can find the stationary points of (AFV) by solving (QFV). 

4) Algorithm for Matrix-Weighted Sum-MSE Minimization: We adopt an alternating minimization 
approach to develop an iterative algorithm for finding the stationary points of (M.FV). In each iteration, 
we focus on determining only one of the sets of parameters Ft, W x , and V while assuming the others 
are fixed. When F x and V are fixed, the optimal linear receive filter at relay k G /C x is exactly w x 1 ^ ISE 
given in ([5]). In addition, as discussed in Section IV-B3 when Ft and W x are fixed, the optimal matrix 
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weights are as follows 

Vf t = V 5fc (E M (F T )),VA:e/Cx. (35) 



It follows from (26) that Efc(Fx,Wx,fc) is a Hermitian and positive semidefinite matrix for any Ft 



and Wx- Combined with (33 1, we have Vgk(E>k,o(FT)) is a Hermitian and positive semidefinite matrix. 
This means that if we always choose = V^ pt according to (35), then is a Hermitian and positive 



semidefinite matrix for k G /Cx- 

What remains is the design of Ft when Wx and Vx are fixed. When Wx and Vx are fixed, we 
need to solve the following optimization problem to determine Ft 



(MFV-F T ): min V tr(V£E fc (F T , W x ,k)) 
f t sf t * — ' 
fce/Cx 

subject to ^ tr(F T , m F T m ) < p Tjfc , Vfc G JC T . (36) 

fj,(m)=k 

Note that (MFV-Ft) is convex with respect to Ft,* for k G /Cx- Let A& > be the Lagrangian multiplier 
associated with the sum transmit power constraint at transmitter k G JCt- Based on the optimality condition 
of {MFV-F^), the globally optimal solution of (MFV-Ft) must has the following form for k G /Cx 

F T , fc (A fc ) = ( H ^(*) Wx, m V m W x , m ff m , Mfc ) + A fc lAT T>fe ] H* Mt) W x ,,V t . (37) 

\mg/Cx / 

Since V k is a Hermitian and positive semidefinite matrix, Eme/C x H m,A 1 (fc) Wx i mVmW x,m H m, At (fc) 
is also Hermitian and positive semidefinite. It can be checked that tr (FT,fc(Afc)F T k (Xk)) is strictly 
decreasing with X k in [0, +oo). The optimal Lagrangian multiplier At > is chosen such that the 
complementary slackness condition of the sum power constraint at transmitter k G /Ct is satisfied. For 
any k G K T , if EmtK x : M (m)=^( F T ) ifc(0)F T ^(0)) < p T ,fc, then F° p * = F T)fe (0). Otherwise, A* is the 
unique solution of the following equation 

£ tr (F T , fc (A fc )F Tifc (A fe )) =p T ,fc. (38) 

mgJCx 
fj,(m)=k 

This equation can be solved by using one-dimensional search techniques, e.g., the bisection method. 

Note that in the proposed algorithm for solving (A4FV), we are able to find the globally optimal solu- 
tions to the corresponding optimization problem in each iteration. Therefore, the algorithm is guaranteed 
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to converge to a stationary point of (MFV), which is also a stationary point of {AFT). Let Ft denote 
the transmit precoders corresponding to the resulting stationary point. It is worthwhile to emphasize that 
it is not guaranteed that we can find a stationary point of (TV). Thus, a two-hop rate mismatch may still 
happen for the resulting suboptimal solution (Ft, Fx) of the original transmit precoder design problem 
(OV). This leaves room for potential improvements in terms of maximizing -R SU m(FT, Fx). 

5) Rate-Matching Transmit Power Control: We propose an iterative power control method for elim- 
inating any residual two-hop rate mismatch corresponding to (Ft, Fx). Let F T k denote the transmit 
precoder for the transmission to relay k G /Cx in iteration n of this method. Note that = F^fc for 
k G /Cx- Let 9^ G K be the norm of F^, i.e., the power allocated for the transmission from transmitter 
fi(k) G /Ct to relay k G /Cx in iteration n. We propose to fix the shapes of the transmit precoders and 
to adjust only their norm 9k for k G /Cx- It follows that f!?1 = 9^ nI T * k n for all n > where ll £ T,k u 

J-A " ||*t,/c||f ||*t,*||f 

has unit norm and represents the shape of Ft,A;- 

We define 6^ = {9^\- ■ ■ , 6*g ) G M. Kx x 1 and 0™ = {9^\ • • • , 9 k % 9^ ■ ■ • , Note that 0^ 
and (9 k ,0_ k ) are used interchangeably in the section. Also, we denote 'H m ,fi(k) = Hm,^(fc) p T, *j| for 
k,m G /Cx- The interference plus noise covariance matrix at relay k G /Cx is rewritten as 

Rx,k(9-l) = ^2 e m^-k,n(m)'Ht,i X ( m ) +°X,k 1 Nx, k , (39) 

while the maximum achievable rate at relay k G /Cx is rewritten as 

fli )fc (* (n) ) = lo g2 det (i + 9 { k n) Ul Ak) [Rx lfc (»5fe)] • (40) 

Thus, it follows that 

= det (i + ^n* kAk) [Rx, k (^-l)} ^Uk,^)) - 1. (41) 
The end-to-end achievable rate corresponding to relay /c G /Cx is written as 

i? fe (0(")) = mm{tR hk (G^), (1 - t) log 2 (l + £ 2)fc )}. (42) 

Note that 1-L* k „fc[Rx,fc(0flfc)] -1, Hfcjifc is independent of 9 k n \ It is also a Hermitian and positive semidef- 
inite matrix. Since det (I + xA) is strictly increasing in x for x > when A is a positive semidefmite 
matrix, both #i, fc (0 (n) ) and £i ife (0 (n) ) are strictly increasing in 9^ when 9 { k ] > 0. 
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Let A^ = {k G K x ■ ti?i, fc (0 (n) ) > (1 - t)log 2 (l + &,*)} be the index set of the relays with a 
dominant first hop in iteration re. If A; € .4.*™) , then excess power is allocated for the transmission to relay 
k. Similarly, let = {k G /C x : tR hk {0 {n) ) < (1 - t) log 2 (l + &,*)} be the index set of the relays 
with a dominant second hop in iteration n. A two-hop rate mismatch happens if and only if A^ ^ and 
£>( n ) ^ 0. When a two-hop rate mismatch happens, we consider an arbitrary kj± G .4.^ and % G ,B( n ). 
It follows ^( n ) n #( n ) = that Au / %. When 0^ is fixed for m G /C x and m / jfe^, since fc^ G „4>), 
there must exist <f>^ G (0,0^) such that 

ti2i, fc (^; 0^)) = (i - 1) log 2 (l + (43) 

Equivalently, it follows that 

det (i + [Rx,fc.(0il)] "^wow) - 1 = tjfcx- (44) 



(n) 

Since the left-hand side of (44 1 is strictly increasing m <p) ', this equation has a unique solution, which 



in) 

can be found by using one dimensional search techniques, e.g., the bisection method. Note that 4> k ^ is 
the rate-matching power for the transmission to relay k A in iteration re. 

The key observation for the proposed power control algorithm is that if a two-hop rate mismatch 

(tI I ll ( 71 I 1 ) 

happens at the end of iteration re, excess power can be reduced by setting 9^. = <j) kA and 0_ k = 
0^l A for an arbitrary k A G A^ . Note that = (1 - *)log 2 (l + = RkJO^), 

Also, this excess power reduction decreases the power of interference observed by all relay m ^ &a, 
leading to i?i im (0 (ri+1) ) > i?i jm (0 (n) ) and hence il m (0 (n+1) ) > R m (6 (n) ). Especially, it follows that 
Ri,k B (d {n+1) ) > Ri,k B ( din) ) and hence Rk B (0 {n+1) ) > R ke {0 (n) ). As a result, the end-to-end sum-rate 
is strictly improved, i.e., Y2keK Rkifl^) < ^2keK x ^fc(^ n+1 ^)- Thus, when a two-hop rate mismatch 
happens, reducing excess power in a controlled manner strictly increases the end-to-end sum-rates. 

Based on the observation, we propose an iterative algorithm for updating the power allocated for the 
transmission from the transmitters to the relays. Specifically, at the end of each iteration re > 0, each 
relay k G /Cr computes £ik(@ ) to check if the transmission to itself is allocated excess power. If 
(,i,k(6^) < i]k, then the power allocated for the transmission to relay k does not need to change, i.e., 
^A™ +1 ^ = @k Otherwise, relay k determines the corresponding rate-matching power cf)^ and feeds back 
the value to its associated transmitter fj,(k) to instruct the transmitter to update 9 k n+1 ^ = 4> k n+1 ^ G (0, 0^). 
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In other words, the power update rule for k G /Cx and n > is 



? (n+l) 



^ n) , if fce^") 5 
otherwise. 



(45) 



This process is repeated until the algorithm converges or the maximum number of iteration is reached. 
It follows from the power update rule that 9^ +1 ^ < for k G /Cx and n > 0. Several properties of 
the proposed algorithm are presented in Remark [5] Remark [6j and Remark [7] 

(n) (n) 

K for /c G /Cx and n > 1. Since 6^ is nonnegative, it is lower 



Remark 5: Note that 



< 



bounded by 0. Thus, it is guaranteed that the proposed power control algorithm converge as the number 
of iterations goes to infinity. 

Remark 6: The resulting solution of the proposed power control algorithm does not causes a two-hop 
rate mismatch. Indeed, let no be the index of the iteration at which the proposed algorithm is convergent, 



i.e., 



for k £ £ x TMs means that £ 1)fe (0M) < Vk for all A; G /c x or A^ = 0. Recall 



that a two-hop rate mismatch happens in iteration n if and only if A^ 1 ' ^ and B^ n > ^ 0. Thus there 
is not any two-hop rate mismatch in iteration uq. 

Remark 7: The proposed algorithm does not decrease R^iO^) over iterations for n > 0. It can be 
showed that if A is positive semidefinite, then det(I + B*(I + xA) _1 B) is strictly decreasing in x for 



x > for any B. From (39) and (41 ), we have 



£i, fc (0 (n+1) ) + l = det 1 + 



> det 1 + 



3 ("+l) 



k,fi(k) 



22 ^k,^{m)^k,tl{m) + a k l 



L me/Cx 
mj^k 



n 



k,fi(k) 



X] ^fe,ft(m)^fe lM (m) + a k l 



m^k 



-1 



n 



k,fj,(k) 



rj k + 1, if k G AW 
(,i,k( ein) ) + !> otherwise, 



(46) 



(47) 



where (|46|) comes from the property that #[, n+1 ^ < 0^ for k G /Cx and (|47|) comes from (|45|) and (|44 1. 



It follows that min{r? fc ,^" fe +1) } > min-T^,^}; equivalently, i4(0 (n+1) ) > R k (9 (n) ) for fc G /C X - 

Remark 8: We can always use the same steps to develop a similar rate-matching transmit power 
control at the relays on the second hop. Nevertheless, because the proposed algorithm is implemented 
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in a distributed manner, it is unclear how to determine when the rate-matching transmit power control 
should be performed on the first hop and when it should be performed on the second hop. 

Note that the proposed power control algorithm does not help find any stationary points of (OV). 
Essentially, it helps eliminate the residual two-hop rate mismatch in (Ft, Fx) to find a potentially better 
suboptimal solution to (OV). 

C. Discussion 

1 ) Summary of the Proposed Algorithm: Recall that the proposed algorithm aims to find high-quality 
suboptimal solutions to the end-to-end sum-rate maximization problem (OV). The development of this 
algorithm is involved with the formulation of a number of optimization problems as illustrated in Fig. 
[2] The flow diagram of the proposed algorithm is presented in Fig. [3j where the notation for the main 
parameters is summarized in Table [j] Note that Phase 1 is a counterpart of Table I in [?]. 

2 ) Distributed Implementation: The proposed transmit precoding algorithm for the relay interference 
broadcast channel allows for distributed implementation. Similar to [?], [?], [?], [?], two assumptions 
are needed. First, each transmitting node has the corresponding local channel state information (CSI). 
Specifically, on the first hop, transmitter k G JCt has the CSI of H m> fc for all m G /Cx; on the second 
hop, relay m G /Cx has the CSI of G qjm for all q G /Cr. Second, there is a feedback channel to send 
information from a receiving node to its serving node, i.e., from receiver q G /Cr to relay x(q) £ ^Qx an d 
from relay m G /Cx to transmitter /j,(m) G /Ct- Specifically, while using the iteratively weighted MMSE 
approach to solve (SV), in each iteration, the receivers need to feedback the updated matrix weights 
and receive filters to the relays [?], [?]. At the end of Phase 1, each receiver m G /Cr computes the 
second-hop received SINR and sends back to its serving relay, i.e. relay x( m )- On the first hop, in each 
iteration of solving (M.FV), the relays also need to send back the updated matrix weights and receive 
filters to the transmitters. 

3) Opportunistic Approach: Recall that it is not guaranteed that the proposed algorithm can find 
optimal solutions to (OV). In fact, the end-to-end sum-rate performance of the resulting solution (Ft, Fx) 
depends on the initial solutions in the first two phases. One method for improving the end-to-end sum- 
rate performance of the algorithm is to use multiple random transmit precoders as initial solutions in 
the first two phases and then to select the best one in terms of end-to-end sum-rate maximization. Such 
an opportunistic approach, however, may require more coordination among the nodes and result in more 
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(d) 



(e) 



» (AFP) — > (GFP) —* (AtFP) ^* (MFT-Ft) 



(f) 



Fig. 2. Illustration of the transformations between different optimization problems to find high-quality suboptimal solutions to 
(OV): (a) the design problem for the transmit precoders at the relays (SV) is obtained from (OV) when the first hop is ignored; 
(b) the design problem for the transmit precoders at the transmitters (FT') is obtained from (OV) and the first-hop performance 
results from the previous step; (c) approximates of the utility function of (FT') given in |24| are used to formulate (AFT); 
(d) (AFT) is equivalently rewritten as (QFT)\ (e) (A4FT) is the matrix-weighted sum-MSE minimization problem that has 
the same stationary points as (QFT)\ the alternating minimization approach is adopted to solve (A4FT) and (f) (MFT-Ft) 
is the design problem to design the transmit precoders at the transmitters in the corresponding iterations. 
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Fig. 3. Flow diagram of the proposed algorithm for the relay interference broadcast channel. 



overhead in the network. 

4 ) Order of Optimization: In the proposed algorithm, the relays are optimized for second-hop sum-rate 
maximization before the transmitters are designed for end-to-end sum-rate maximization. Note that the 
rates on two hops are matched with each other per relay in the design of the transmitters. In particular, 



when the allocated first-hop rates for k G /Cx and m S /Cr are determined as in (15 1, at each 
relay k, the normalized first-hop rate i-Ri ^F-p) is matched with the sum of normalized second-hop rates 
(1 — i)-R2,fc,sum(Fx)- There is another order of optimization where the transmitters are optimized for 
first-hop sum-rate maximization before the relays are optimized for end-to-end sum-rate maximization. 
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TABLE I 

Summary of notation 



Notation 


Parameters 


X(k) 


index of the relay aiding receiver G /Cr 




index of the transmitter aided by relay k G /Cx 




precoder at relay G /Cx for transmission to receiver q G /Cr 


F X ,g 


resulting transmit precoder at relay x(q) G /Cx for transmission to receiver q G /Cr 




receive filter at receiver q G /Cr 


u g 


matrix weight for MSE at receiver q G /Cr in the second-hop weighted sum-MSEs 




effective SINR corresponding to second-hop sum-rates from relay k G /Cx 


m 


rate-matching SINR at relay k G /Cx on the first hop 




precoder at transmitter n{k) G /Ct for transmission to relay k G /Cx 


Ft,A: 


resulting precoder at transmitter /x(/c) G /Ct for transmission to relay k G /Cx in 
the second phase 




receive filter at relay k G /Cx 


v fc 


matrix weight for MSE at relay k G /Cx in the first-hop weighted sum-MSEs 




Frobenius norm of Ft,£; in the third phase, i.e., the power for transmission from 
transmitter /j,(k) G /Ct to relay k G /Cx 




effective received SINR at relay k G /Cx 



In this case, the second-hop achievable rates i?2,m(Fx) for m G /Cr are determined given the knowledge 
of the first-hop achievable rates Ri^(Ft) for k G /Cx- For a fixed first-hop rate allocation strategy (3k, m 



for k G /Cx and m G /Cr, the end-to-end sum-rates in (111 can be rewritten as 



#2,fc,sum(F x ) = ^ min{/3 x(m)im ,(l-t)i? 2 ,m(Fx)}. (48) 

mg^x 

Note that the end-to-end sum-rate performance in this case depends significantly on the selected first-hop 
rate allocation strategy. Moreover, when the set of (3k,m for k G /Cx and m G /Cr is given, the same steps 



as in Section IV-B can be used in the design of relays in this case if matching the rates on two hops 
is performed per receiver. What remains is to determine the optimal first-hop rate allocation strategy in 
terms of end-to-end sum-rate maximization. Unfortunately, it is challenging to do this before designing 
the relays, i.e., without the knowledge of second-hop achievable rates i?2,m(Fx) for m G /Cr. This is 
the main reason that we decided that the relays are designed before the transmitters. In the special case 
where each relay serves a single receiver, the two orders of optimization are the same since there is no 
need for first-hop rate allocation. 
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V. Simulations 

This section presents Monte Carlo simulations to investigate the end-to-end sum-rate performance of 
the proposed algorithm. We consider only symmetric relay systems with k = -Wt and pr,fc = Pi for 
k G /C T ; Nx >m = N x , d ljTn = d\, and px, m = Px for m G /C x ; and iVR i(? = Ar and d 2>q = d 2 for q G /Cr. 
Each transmitter is aided by Kx/Kt relays. Each relay forwards data from its associated transmitter to 
Kr/Kx receivers. We denote the system as (iV^ x N-£ x X N^ B ,d\ x d 2 ). Except when it is stated 
explicitly, we use t = 0.5. The channels are flat both in time and in frequency. The channel coefficients 
on two hops are generated as i.i.d. zero-mean unit-variance complex Gaussian random variables. Path 
loss is not considered in the simulations, thus the average power values of all cross-links on the same hop 
are equal to each other. The power values are normalized such that p^ — px> &x. m — 

1 for m G /Cx and 

<7R,g = 1 for g G /Cr. The plots are produced by averaging over 1000 random channel realizations. In 
each channel realization, the initial transceivers are chosen randomly. The maximum number of iterations 
in the first two phases is 2000 while that of the last phase is 30. 

A. Convergence of the Proposed Algorithm 

Note that if each of its three phases is convergent, the proposed algorithm is also convergent. The 
first two phases are based on the prior work for the single-hop MIMO broadcast channel, of which the 
convergence is well validated in [?]. Fig. [4] shows the convergence behavior of the proposed rate-matching 
transmit power control method for a channel realization of the system (2 3 x 2 6 x 2 12 , 1 x 1). We observe 
that this method converges in few iterations and it does so monotonically. 

B. Benefits of Each of the Last Two Phases 

In this section, we investigate the benefits of each of the last two phases by comparing the proposed 
algorithm with three alternative algorithms. Fig. [5] shows the block diagrams of all four algorithms. All 
algorithms use the same design of the second-hop precoders, which is the first phase of the proposed 
algorithm. Although they use different approaches to designing the first-hop precoders, they always use 
the same initial solution. The naive approach is used as the main baseline strategy, which is labeled as 
Baseline. The output of the proposed algorithm is labeled as Final output. The output of the second 
phase in the proposed algorithm is labeled as After phase 2. Since the proposed rate-matching power 
control method for first-hop precoder adjustment can be applied for any (Ft, Fx), we proposed an 
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Fig. 4. Convergence behavior of the rate-matching power control algorithm for the (2 x 2 6 x 2 12 ,1 x 1) system. 



algorithm that applies the proposed rate-matching power control for first-hop precoder adjustment to the 
output of Baseline. We label it as Baseline + PC. Note that the comparison between After phase 2 and 
Baseline gives us insights into the benefits of the second phase of the proposed algorithm. Similarly, the 
comparison between Final output and After phase 2 and that between Baseline + PC and Baseline give 
us insights into the benefits of the third phase of the proposed algorithm. 
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Fig. 5. Block diagrams of the proposed algorithm and the three reference algorithms. The dashed lines represent the links used 
for exchanging design parameters. The solid lines represent the output of the corresponding blocks. 



Fig.[6]shows the average end-to-end sum-rates achieved by the four algorithms for the (2 3 x 2 6 x 2 12 , 1 x 
1) system. All algorithms achieve non-zero end-to-end multiplexing gains thanks to their capabilities of 
interference management on two hops. The proposed algorithm outperforms the other algorithms for all 
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Fig. 6. Investigation of the benefits of the last two phases for the (2 x 2 6 x 2 12 ,1 x 1) system. 

transmit power values. At pt = Px = 30 dB, Final output provides a gain of 10% over After phase 2, 
which in turn provides a gain of 44% over Baseline. Also, Baseline + PC provides a gain of 28% over 
Baseline. The results mean that thanks to the consideration of t and £,2,k for k G Kx in the design of Ft, 
the second phase of the proposed algorithm is able to alleviate a two-hop rate mismatch to obtain higher 
end-to-end sum-rates. In addition, the results show that there exist two-hop rate mismatches at the output 
of both After phase 2 and Baseline. This is expected since the second phase of the proposed algorithm 
aims only at maximizing an approximate of the end-to-end sum-rates. We notice that the two-hop rate- 
mismatch of the output of Baseline is relatively more significant than that of the output of After phase 2. 
This emphasizes the benefits of two-hop rate-mismatch alleviation of the second phase of the proposed 
algorithm. Note that After phase 2 outperforms Baseline + PC in this simulation scenario, however, it is 
not guaranteed that After phase 2 always outperforms Baseline + PC. 

Note that the relative gains provided by each of the last two phases depend on the configurations. 
Fig. [7] provides the simulation results for After phase 2 and Baseline for the following three systems: 
(2 3 x 2 6 x l 12 , 1 x 1), (2 4 x 2 4 x 2 4 , 1 x 1), and (4 4 x 4 8 x 2 16 , 2 x 2). We observe that even before 
rate-matching power control, the proposed algorithm always outperforms the baseline. At pt = Px = 30 
dB, the gain is 60% for the system (2 3 x 2 6 x l 12 , 1 x 1), 22% for the system (2 4 x 2 4 x 2 4 , 1 x 1) and 
40% for the system (4 4 x 4 8 x 2 16 , 2 x 2). 
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Fig. 7. Comparison of average end-to-end sum-rates of After phase 2 and Baseline to show the benefits of the second phase 
of the proposed algorithm. 

C. Opportunistic Solutions 

We adopt an opportunistic approach to improve the end-to-end sum-rate performance. Let N be the 
number of random initializations in the opportunistic approach. The proposed algorithm is repeated 
iV times with the random initializations and choose the one with the highest achievable end-to-end 
sum-rates. Fig. [8] provides the average achievable end-to-end sum-rates of the opportunistic solutions of 
the proposed three-phase algorithm with the number of initializations N G {1,2,5,25} for the system 
(2 3 x 2 6 x 2 12 , 1 x 1). To provide a benchmark, we also show the results for the fixed average normalized 
second-hop sum-rates, which provides an upper-bound for the solutions. As expected, increasing TV 
improves the end-to-end sum-rates achieved by the proposed algorithm with opportunistic implementation. 
At Pt = Pr = 30 dB, the opportunistic solution with N = 5 nearly doubles the end-to-end sum-rates 
when compared to the baseline. Nevertheless, the additional gains obtained by using an additional random 
initialization in the opportunistic solutions decreases in N. At Pt = Pr = 30 dB, the opportunistic 
solution with N = 25 achieves nearly 75% the value of the upper-bound. 

D. Varying Timesharing Values 

Fig. [9] presents the average achievable end-to-end sum-rates as functions of the timesharing values for 
the system (2 4 x 2 s x 2 16 , lxl) for the following two cases: i) Case 1 with pt = px = 30 dB and ii) 
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Fig. 8. Average end-to-end sum-rates of opportunistic solutions of the proposed three-phase algorithm for N G {1,2,5,25} 
for (2 3 x 2 6 x 2 12 , lxl). 

Case 2 with = 30 dB and px = 20 dB. We notice that for the proposed algorithm and the baseline in 
each case has the same optimal timesharing value, t\ = 0.5 for Case 1 and t2 = 0.425 for Case 2. These 
optimal timesharing values approximately equalize the average normalized sum-rates on two hops. This 
emphasizes the importance of matching the rates on two hops. In addition, we observe that thanks to the 
explicit consideration of t for matching the rates on two hops, the proposed algorithm has large gains, 
between 50% and 70%, over the baseline for t e [0.1,0.9]. 

VI. Conclusions and Future Work 

We proposed a cooperative algorithm for jointly designing transmit precoders at the transmitters and 
relays for the decode-and-forward relay interference broadcast channel to (approximately) maximize 
end-to-end sum-rates. The naive approach of applying single-hop interference management strategies 
separately for two hops of the system causes rate mismatch, leading to low end-to-end sum-rates. The 
main challenges were interference management and rate mismatch between the rates on the two hops. Our 
solution is a three-phase cooperative algorithm. The first phase focuses on the design of precoders at the 
relays to maximize the second-hop sum-rates. The second phases uses the knowledge of the timesharing 
value and the second-hop achievable rates to design the precoders at the transmitters to maximize an 
approximate of the end-to-end sum-rates. The last phase uses power control at the transmitters to eliminate 
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any residual rate mismatch to further improve the end-to-end sum-rates. The proposed algorithm allows 
for distributed implementation and has fast convergence behavior, making it suitable for practical systems. 
Simulations showed that the proposed algorithm obtains much higher end-to-end sum-rates than the naive 
approach. 

The work assumes instantaneous and perfect CSI is available at the transmitters and relays. Of course, 
perfect CSI is not realizable in practical systems where there is channel estimation error, time-variation 
in the channel, thermal noise, errors on the feedback link, and signal processing delay. Future work could 
investigate and quantify impacts of imperfect CSI, e.g., due to channel estimation or CSI feedback delay, 
on the proposed algorithm. Based on this knowledge, future work should focus on developing algorithms 
that are robust in presence of CSI uncertainty. Our proposed algorithm provides benchmarks for future 
work that accounts for more practical impairments. 



